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Abstract. We consider the total cost of cutting down a random rooted tree 
chosen from a family of so-called very simple trees (which include ordered trees, 
d-ary trees, and Cayley trees); these form a subfamily of simply generated 
trees. At each stage of the process an edge is chose at random from the tree 
and cut, separating the tree into two components. In the one-sided variant 
of the process the component not containing the root is discarded, whereas in 
the two-sided variant both components are kept. The process ends when no 
edges remain for cutting. The cost of cutting an edge from a tree of size n 
is assumed to be n". Using singularity analysis and the method of moments, 
we derive the limiting distribution of the total cost accrued in both variants 
of this process. A salient feature of the limiting distributions obtained (after 
normalizing in a family-specific manner) is that they only depend on a. 



1. Introduction 

Consider the following process on a rooted tree with n vertices. Pick an edge 
uniformly at random and "cut" it, separating the tree into a pair of rooted trees; 
the tree containing the root of the original tree retains its root while the tree not 
containing the root of the original tree is rooted at the vertex adjacent to the edge 
that was cut. In the one-sided variant of the problem the tree not containing the 
original root is discarded and the process is continued recursively until the original 
root is isolated. In the two-sided variant the process is continued recursively on 
each of the rooted trees. Assume that the cost incurred for selecting an edge and 
splitting the tree is tn- In this paper we derive the limiting distribution of the total 
cost accrued when the tree is a random very simple tree (defined below) and tn = n" 
for fixed a > 0, for both the two-sided variant (Theorems 14 . 41 \TM and IT77|l and the 
one-sided variant (Theorem I5.1|l . A salient feature of the limiting distributions 
obtained (after normalizing in a family-specific manner) is that they only depend 
on a. 

In the one-sided variant, the case i„ = 1 (i.e., a — 0) corresponds to the num- 
ber of cuts required to disconnect the tree. For this random variable, Meir and 
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Moon ^5 derived the mean and variance for Cay ley trees; Chassaing and Mar- 
chand derived the limiting distribution for Cayley trees. Panholzer obtained lim- 
iting distributions for non-crossing trees |16| and very simple families of trees jl5j . 
Recently Janson extended these results to all simply generated families |10| . 

The interest in the two-sided variant stems from the fact that when the very 
simple family is Cayley trees, the process is equivalent to a probabilistic model 
(the "random spanning tree model") involved in the Union-Find (or equivalence- 
finding) algorithm. Knuth and Schonhage derived the expected value of the 
cost in the cases (among others) i„ ~ a^pn and t„ = n/2. These results were later 
extended ^3 to the cases i„ = when a > 1/2 and t„ — 0{n") when a < 1/2. 
(Some of these expected values were rederived using singularity analysis in jj.) 
In (2| , Chassaing and Marchand derive limit laws for the costs considered by Knuth 
and Schonhage. 

We treat both variants of the destruction process using singularity analysis [5], 
a complex-analytic technique that relates asymptotics of sequences to singularities 
of their generating functions. We rely on applicability of singularity analysis to the 
Hadamard product (the term-by-term product) of sequences 0] and the amenability 
of the generalized polylogarithm to singularity analysis T . 

The organization of the paper is as follows. In Section |21 we define families 
of very simple trees, noting the key "randomness-preservation" property that is 
crucial for the application of our methods. Section|21establishes notation and other 
preliminaries that will be used in the subsequent proofs. In Section0]the two-sided 
variant is considered, and Section ideals with the one-sided variant. 

Notation. In the sequel we will use In to denote natural logarithms and log when 
the base of the logarithm does not matter. 

2. Very simple trees 

An ordered tree is a rooted tree in which the order of the subtrees of each given 
node is relevant. Given a sequence (0i)i>o of nonnegative numbers (called a degree 
generating sequence) with (/iq = 1, a simply generated family T of trees is obtained 
by assigning each ordered tree T the weight 

where d(v) is the outdegree of the node v. Let JF„ denote the set of trees in T 
with n nodes, and let T„ denote the weighted number of trees in JF„, i.e., 

A random simply generated tree of size n is obtained by assigning probabil- 
ity w(T)/Tn to the tree T G Tn. Many combinatorially interesting families such as 
(unweighted) ordered trees, Cayley trees, Motzkin trees, and d-ary trees are simply 
generated. It is also well known that simply generated trees correspond to certain 
conditioned Galton-Watson trees; see the introductory section of ^Hl for the pre- 
cise connection. It is well-known that the generating function T(z) := X]ti>i-^"'^" 
satisfies the functional equation 

T(z) = z$(T(z)), 
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where := J2k>o't'k't'' is the degree generating function of the family. For 

further background on simply generated trees we refer the reader to |14j . 

In this paper we consider the subclass of simply generated families, called very 
simple families, that, among simply generated families, are characterized by the 
following property. 

Choose a random simply generated tree from the family J^n and then one 
of its n — 1 edges uniformly at random. Cutting this edge produces a pair 
of trees of size k (the one that contains the root) and n — fc, as described 
in Section ^ Then the subtrees themselves are random simply generated 
trees from the family J^k and Tn-k- 

It is clear that the "randomness-preservation" property of very simple trees al- 
lows a simple recursive formulation [see H3.1|l and (|3.12l) ] of the total cost of de- 
stroying such a tree. 

Panholzer |15l Lemma 1] characterized the degree generating functions of very 
simple trees; the relevant constraints are summarized in Tabled 

2.1. Singular expansions. As is usual for treatment of simply generated families, 
let T denote the unique root of t^'{t) = with < t < R, where R is the radius 
of convergence of the series $. Let p := r/$(T). Let Z :— l—p~^z, and let A denote 
a generic power series in Z , possibly different at each occurrence. Then as z —^ p, 
the dominant singularity for T[z), a singular expansion for T is 9, Theorem VII. 2] 



(2.1) T{z) - T - bp^/^Z^/'^ + ZA+ Z^/^A, 



where b := $(r)y -^r^^^yj- immediately from singularity analysis that 

(2.2) Tn ^ cp-"n-3/2(i + n-^U), 

where c :— 6^^/2/(2^/7?) = [$(r)/27r<I>"(T)]^/2 the sequel we will also use 

(2-3) .= r 

Differentiating the expansion (|2.1|l tcrm-by-term Ti, Theorem 6] we get 
T\z) ^ + ^ + z^'^A. 



Since z = p — pZ, 



zT'(z)^^p'^'Z-'^'+A+Z'^'A. 



(2.4) zT'iz) ~ -p 

The constants ap and ai described by Table ^ are fundamental constants for our 
analysis; see especially (13.211 . Using (12. If) and H2.4|) we get 

(2.5) 1 + 2aoT{z) + aizT'{z) - a^p^/^^Z-^/^ + A+ Z^/^A 
and 

1 - aiT[z) - (1 - air) + aibp^'"^ Z'^''^ + ZA+ Z^/'^A. 
It is easily verified that for each very simple family 1 — aiT = (this fact will be used 
numerous times in subsequent calculations), so that the constant term vanishes in 
the singular expansion of 1 — aiT{z). This leads to 

z[l - aiT(z)] - j^zA+ Z'^I'^A 
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and consequently 

(2.6) z-i[l - aiT(z)]-i ^ ^ _^ ^ ^1/2^^ 
We will also need 

(2.7) + + 

1 yz) 

which follows from (|2.1fl . 



3. Preliminaries 

Throughout, — denotes equality in law (or distribution) and — > denotes con- 
vergence in law. Recall that the Hadamard product of two power series / and g, 
denoted by / © 5, is the power series defined by 

(/ © g){z) EE f{z) © giz) ^ /„ff„z", 

n 

where 

f{z) = fnz" and g{z) = ^ g^z". 

n n 

3.1. Two-sided destruction. The cost of cutting down a very simple tree of 
size n, call it Xn, satisfies the distributional recurrence 

(3.1) X„ = Xk,^ + X*_j^^^ + 1„, n > 2; Xi=ti, 

where tn, for n > 2, is the toll for cutting an edge from a tree of size n. Here Kn, 
the (random) size of the tree containing the root, is independent of {Xj)j>i and 
(X*)j>i, which are independent copies of each other. The splitting probabilities are 
given by 

(3.2) Pr[i^n = fc] Pn,fc = (oifc + 00)7^^— fc = l,...,n-l. 

[n - 1)T„ 

TableHgives the constants oi and qq for each type of very simple family; see (14)- 
(16) in IIS]. Here at := ^i+i/c/ii, i = 0,1, where {4>i)i>o is the degree generating 
sequence of the simply generated tree. It is easy to check that family A is Cayley 



Family Generating function Constraints ai oq 
A ^ ^0 

C [1 — {2ai — Q;o)t]~ 2°i-°o 2q;i — OfQ > 2ai — (2q!i— ao) 

Table 1. Generating functions for very simple families. For each 
family, ao > is also a constraint. 



trees, family B is d-ary trees, and family C contains unweighted ordered trees. (As 
it turns out, the distributional recurrence for Cayley trees is identical to the one 
obtained for the Union-Find process studied in ^21E1E] — see Remark ^31 below.) 
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Define :— EX^. Taking sth powers of both sides of H3.1|l and taking expec- 
tations by conditioning on Kn, we get 



n-l 



(3.3) 

where 

(3.4) 



fc=i 



E 

Sl+S2+S3=S 



Si, S2, S3 



[S2],,[ss] 

'n—k ' 



[si 

and ii\ — ti- Define generating functions 

n>l ""-^l n>l 

[Observe that /x[°l(z) = T(z).] Multiply by (n - l)T„z" and sum over n > 2. 
The resulting left side is 

J2{n - l)Tn^^^^^z- = E(" - = ^5.A^W(z) - A'W(^), 

n>2 n>l 

where 9^ denotes derivative with respect to z. Similarly, the resulting first term on 
the right side is 

ai [2 (9.m'''(2)) T{z) + zr'(z)A*W(z)] + 2aoH^'\z)T{z). 
The resulting second term on the right side is 
rW(z) := E(" - = E(" - ^)Tnr^:^^" 



(3.5) 



E 



S2.S3<S 



Si, S2, S3 



i®^^ (z) [aiz (z)) Ail^^'l (z) + aoM^''' (^)m'''' (^) 



Thus (|3.3() translates to 

za^^W(z) - ^W(z) 

= ai 

i.e., 

(3.6) 
where 

(3.7) 
and 
(3.8) 



z (a.AiW(^)) Tiz) + zr'(z)AiW(z)J + 2ao/^M(z)T(z) + rW(z), 

9.Mt^l(z)+p(z)A^W(^)=ff'''(^), 
1 + 2aoT(z) +aizT'(z) 



p(z) - 



5W(z) := 



z[l — air(z)] 

rM(z) 
z[l-aiT(z)]' 
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with /^t"! (0) = 0. By variation of parameters (see, for example, Jl , 2.1-(22) and Prob- 
lem 2.1.21], the general solution to the first-order linear differential equation H3.6|l 
is given by 



(3.9) 

where 

(3.10) 



exp 



p{t) dt 



^W(z):= / gW(t) 



I exp 



p{u) du 



dt + P„ 



with zo chosen as follows and f3s an arbitrary constant. 

The integrand p{z) defined at H3.7|l and appearing in (|3.9|l - (|3.1()(l is asymp- 
totic to —1/z as z — > and has [see (|4.1|) below] another singularity at z = p. 
In H3.9|) ^ (|3.1U|I we may choose (and fix) zq arbitrarily from the punctured disc of 
radius p centered at the origin. Then, in H3.10|l , as t ^ we have 



exp 



p{u) dt 



exp 



exp 



I du 

zo \ 

— In < + In Zq 



p{u) + - 



p{u) + - 



ZQe"'t ^, where a :— 



r 



p{u) + - 



1 

u 

1 

u 

1 

u 



du 
du 
du 



whereas, using 1)3. 8() and 1)3. 5(1 . 



rW(t) 



thus the integrand in (|3.10() has no singularity at t = 0. 

Now we obtain the particular solution of interest, using the boundary condition 
/il^l (z) ~ t\Tiz as z — > 0. We find the constant j3s is specified in terms of zq as 

(3.11) Ps ^ zoe''t{T^. 

Remark 3.1. One can check for each very simple family that 

1 + apr 

and for any simply generated family that 



T'{z) = 



1 - 



<i"(T(z))T(z) 



Thus 



p{z) 



T(z)[l-aiT(z)] 



1 + 2aoT{z) + 



aiT{z) 



1 - 



(aQ+ai)T{z) 
l+aoT{z) 
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This leads to 
p{t) dt = - 



1 



T{z) 
T{zo) "^(1 - 



T{z) 



T{zo) 



Go 



1 + 2aoT 



ai 



aiT 



{aa+ai)T 
l+aoT 



(ap + ai)T 
l + aoT 



dT 



1 



T 1 + aoT 1 — aiT 



dT 



In 



1 — aiT[z) 



T(z)(l + aoT(z)) 



-In 



1 - aiT(zo) 



r(zo)(l + aor(zo)) 



and finally, again using the boundary conditions on /^I'^l (z) as z — s- 0, to the following 
explicit form of (|3.9|) : 



,] , _ r(z)[l + aoT(z)] 



1 - aiT{z) 



9" i ^T7 , dt + tf 



3.2. One-sided destruction. Here, the cost of cutting down a very simple tree 
of size n, call it Ym satisfies the distributional recurrence 

(3.12) r„^yK„+t„, n>2- Yi=h, 

where for n > 2, is the toll for cutting an edge from a tree of size n and the 
splitting probabilities are given by at H3.2|l . 

Defining := EF^, one obtains from equation (|3.12(l by conditioning on Kn 
the recurrence relation 

n-l 

(3.13) =5]p„,fc4^1+rW, n>2, 



k=l 



where 



Sl+S2=S, 
S2<S 



n-l 



S2] 



k=l 



[si . ... 

and = tf. Using the same notation as in Section l3. II we obtain the following 
differential equation by multiplying (|3.13l) by (n— l)T„z" and summing over n > 2: 

z9,/iW(z) - A^W(z) = T{z){alz^z^l^'\z) + ao^i^'\z)) + r^'^z), 

where 

(3.14) rW(z):= 

4 

s 

This can be written as 



S2<S 



Si 



t®^! (z) [T(z) (aiza./il^^l (z) + ao/.[^^l (z))] . 



(3.15) 
with 

(3.16) 



a.A^W(z)+p(z)A.W(z)=gW(z), 



l + aoT(z) [,] rM(z) 

z[l — ail [z)\ z[l — ail [z)\ 



One can check that for each very simple family, p{z) = —dz ln(T(z)), so that we 
obtain as general solution of the first order linear differential equation H3.15|l : 

MW(z)=T(z)^'^dt + Cr(z), 
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and finally by adapting to the initial condition dzfJ-^'^^{z)\^_^ ~ Tijil = Titf, that 
the integration constant is given as C = tf . Therefore, we get 



(3.17) 



^l^\z) = T{z)£^dt + tlTiz). 
4. Two-sided destruction 



We begin by obtaining a singular expansion for p{z) at H3.7|) . Using (|2.5|) 
and in (|X7jl we get 



(4.1) 

Integrating this singular expansion term-by-term 01 Theorem 7], 

p{t) dt 



'^InZ-^ +A+ Z^^^A. 



Thus 

(4.2) exp - f p{t)dt 
where 

(4.3) e:=(l-p-izo)'/'exp 
Taking the reciprocal of 1)4.2(1 gives 



(4.4) 



exp 



p{t) dt 



'^Z-^'^ + A + Z^I^A, 



p{t) + '—{i-p-h)-^ 



C^Z^/'^ + ZA + Z^/^A. 



dt 



Let us now consider two-sided destruction with the toll t„ — n", with a > 0. 
(Notice that the case a = is trivial since then the total cost of destruction is simply 
the number of edges in the tree, which is always n—1.) The toll generating function 
t[z) is the generalized polylogarithm \A-a,o{z), which is amenable to singularity 
analysis 7, Theorem 1]. 

4.1. Expectation. Now we obtain a singular expansion for r[^l(z) defined at H3.5|l . 
recalling that ^["l(z) = T{z): 

(4.5) r[il(z) = t{z) [aizT'{z)T{z) -j- a^T^{z)]. 

Using ((2.1|l we conclude that 

T^(z)^ A + Z^/'^A, 

and using H2.4|) that 



(4.6) 



a^zT'{z)T{z) + aoT\z) ^ + _4 + z^/2j^_ 



We will use the Zigzag algorithm of ^ to obtain a singular expansion for r^^\z). 
We recall the use of the notation A/" to denote a generic power series in 1/n, possibly 
different at each occurrence. By singularity analysis, 

(4.7) p^^[z-][a^zT'{z)T{z) + a^T\z)] ^ p''^^"^ + n-^'^N. 

2 J-K 
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Thus 

(4.8) p"[z"]rW(z) ^ + n"-iAA. 

2 V'"' 

Until further notice, suppose a ^ {i, |, . . .}. Then a compatible singular expan- 
sion for rl^l (z) at (|4.5|l is obtained as 

(4.9) rW (z) ~ pi/25£(^±i) ^-a-i + + A 
Recalling (|3.8|) and H2.6|) we have 

(4.10) <?W(z)- + + + A 

p 2V7r 

Using this expansion and (|4.4|) . 



g[^l(z)exp 



p(i) 

2Q 



20^ 

By (IH-lOf) and Theorem 7 of we may integrate this expansion term-by-term to 
get a complete singular expansion for A. IfQ!^{l,2,...}, we have 

C 2V7r 

where is a constant. [The value of Lq is immaterial unless < a < 1/2, in which 
case sec (14. 14(1 .] On the other hand, if a G {1,2,...}, a logarithmic term appears 
upon integration, so that 

AW(z) ~ - ~ + Z-°+U + + In 

? 2V7r 

where is a constant. Combining these expansions with H4.2|l . we finally obtain 
[recalhng 

(4.11) (z) ^ T ^^"~^^ Z~° + + ^—+1^ + + + 

2v7r 

when a ^ {1, 2, . . . } and 
(4.12) 

/xW(z) - 2) ^-0 ^ ^-a+i^ + + Z-i/2(lnZ-i)^ + (In Z-i)^ 

2V7r 

when a G {1, 2, . . . }. Note that the remainder in 14.12|) is 0(|Z|^"+5) unless a = 1, 
in which case it is 0{\Z\-^ InZ^^) = 0{\Z\-^~^) for any e > 0. 
When a > 1/2 and a^{l,2,...}, by singularity analysis we have 

2V7rr(a) 

so that, recalling (|2.2() and (|2.3|) . 

v2r(a) 
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When a £ {1,2,...}, starting from H4.12|l and the note fonowing that display, we 
can similarly derive the expansion 



-n°+2 +0(n") + 0(n log n). 



\/2r(a) 

When < a < 1/2, a similar computation yields 



(4.13) 
where 
(4.14) 



V2r(a) 



Lq I g'^l(t)exp I p{u)du 

Jo Uzn 



dt + l3i, 



with p and g^^' defined at (|3.7|) and H3.8|l . respectively, and ^ and /3i at (|4.3|) 
and H3.11|l . respectively. 

When Q!€{|,|,...}, one can check that logarithmic terms appear in the singu- 
lar expansion compatible with (|4.8|l but the lead-order term and asymptotic order 
of the remainder are unchanged. Indeed, now 



(4.15) 

and consequently 



^W(z)=r^^i^^Z-" + 0(|Zr"+^) 



V2r{a) 

Finally we consider a = 1/2. Now, a compatible singular expansion for H4.8|l is 



l(z) 



2^ 



Z-^ + {log Z)A + A. 



Proceeding as in the case a ^ 1 /2 we have the singular expansions 



(z) ^ — + z-^A+ (Z-i/2 log Z)A + Z-^I^A + (log Z)A, 

2/3 



5W(z) exp 



V{t) dt 



^[11 (z) - In + Li + + + log Z)A + (z3/2 log Z)A, 

4 2^71 



r 1 



Z-^ +Z-^A+ (log Z)yi + + (Zi/2 log Z)yl, 



where 



(4.16) ii ^"13^(0 exp 



p(u) 



T 1 



(l-p^h) } dt + Pi. 



This leads to 
r 



mW(z) 



2V^ 



In ^-1 + ^LiZ-^l^ + (log Z)yl + + (Zi/2 log Z)A + 



so that by singularity analysis and (|2.2(l we have 



(4.17) V2^ 
where a is defined at 1)2. 



nlnn 



e^i 



CVTT 



(7 + 21n2) 



+ (?^i/2 log 7i)AA + (log n)N + J\f, 
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4.2. Higher moments and limiting distributions. We proceed to higher mo- 
ments. We wih consider separately the cases a > 1/2, a < 1/2, and a = 1/2. 
We present the details for a > 1/2 and sketch the main ideas for the other cases. 
Throughout a' := a + 

Proposition 4.1. Let a > 1/2 and e > 0. Then 



whe 



with 



min{a — ^,^} ifa^^l 
i — e if a = 1 



Ci 



20F ' 



and, for s > 2, 
(4.18) 

C, = p-V25-i 



1 

/ - 1 ^ 

k=l 



feOs_fe + ST 



r{sa' - 1 



r((s-iK-i) 



Proof. The proof is by induction on s. The claim is true for s = 1 by H4.11|l . 14.12|1 . 
and H4.15|l . Suppose s>2. We analyze each term in the sum for r[*l(z) at H3.5I) . 
If both S2 and S3 are nonzero, then by the induction hypothesis, 

zd^fi^'^^ (z) = Cs, (^S2a' - + 0(1^1-^^"'-^"), 

so that 

Also, = 0(|Z|-(''2+''^')"'+i). jjgj^pg 

aiz (a^Ai'"'' (z)) A^'"'! (2) + aoA*'"'' (^)m'''' (-2) 

Taking the Hadamard product of this expansion with t®''^ (z) (using the Zigzag 
algorithm again) gives the contribution of such terms to A^^ (z) as 



Sl, S2, S3 



aiCsjCs, S2a - - 



, i\ r(sa'-f) 



2; r((s2 + S3)a') 



Notice that if si 7^ the contribution is 0(|Z|^*" +^). 

Next consider the case when S2 is nonzero but S3 = 0. By the induction hypoth- 
esis and the singular expansion of T at H2.1|l . 

z (a,M''^l(^)) T(z) = tCs, (^S2a' - Z-^^"'-^ + OdZI-^^^'-H-?). 
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Also /i["=^l(z)T(z) ^ 0{\Z\-^^°-'^^). Hence 
aiz (9,//^'^' (^)) T{z) + aoAi'^^l {z)T{z) 



Taking the Hadamard product of this singular expansion with t®^^ (z) we get that 
the contribution to r[*l(z) from such terms is 

sij \ 2y r(s2a' + i) 

Notice that si > 1 and that when si > 1 the contribution of such terms is 
0(|Z|-^"'+^). 

We move on to the case when S2 = but S3 is nonzero. By the induction hypoth- 
esis, (|2IU> and ITU), we have T{z)p\^'^\z) = 0(|Z|-^3"'+5) and zT' {z)^l^''^\z) = 
0(|Z|-''3"'). Thus 

aizT'{z)y.^'^\z) + aoT{z)^l^''\z) = 0(|Zr^^"'). 

Taking the Hadamard product with t®'^i(z) we see (recalling si > 1) that the 
contribution to r['*l(z) from these terms is 0(|Z|^'*" +2). 

Finally we consider the case when S2 — — 0. In this case, using H4.6|l it is easy 
to verify that the contribution to r[*l(z) from this term is 

Summing all the contributions we see that 



(4.19) 
where 



Ds := ai 



(z) = DsZ-'°'' + 0(|Z|-""'+«), 



.k=l 



r((s-l)a' + i) 



Thus, using and ifTH]! . 
whence, using (|4.4|l . 

= rV^^^^a^^^-^D^Z-""' + 0{\Z\- 



5'''<(2;)exp 



To get A^^^ (z) at H3.10|l we integrate this singular expansion, noting that since s > 2 
and a' > 1, we have sa' > 2. Hence 

sa' — 1 

Now by (133 and g^J, 

^W(z) = p-i/2„-ij,-i_^^— '+1 + 0(|Z|-^«'+H9). 
sa' — 1 
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Taking 



sa' — 1 



P 



sa' — 1 



s-l 

E 

Lfe=l 



ka' - -] CkCs-k + ST 



[(s-l)a'-i]r(sa') 



r((s-l)a' + i) 



7— T^'^«-i 



r(sa' - 1) 



r((s-l)a'-i) 



completes the proof. □ 

[si 

Using singularity analysis we can now derive asymptotics for the moments fin ■ 
Theorem 4.2. Let a > 1/2. Then, as n — > oo, 



where :— o"'^ "t-s (which does not depend on the very simple family) is 

given by 

r(a-i) 



mi 



V2T{c 



and, for s > 2, 
(4.20) 

s-l 



rus 



l_ fs\ V{ka' -\)V{{s-k)a' -\) 



fc=i 



T{sa' - i) 



.sr(sa' - 1) 
rukms-k H — : t-^Jt-s-i- 



Proof. Using singularity analysis and ProDOsition l4.1l 



77 2 

pVWt„ = c , + OK" -5-'?), 

1 (sa' — 2) 



and using the asymptotics of T„ at H2.2|l . 



C, 



cr{sa' - i) 



Then 



where 
(4.21) 

Thus, using 2y/nca = V^t, 



C, 



mi 



cV{sa' 



Ci r(a - i) 



caT{a) ^/2T{a) 
Using (|4.18(l . 14.21|l . and the identities 
1 



cp 



-1/25-1^^ = 1 ^ r(2; + l) = xr(a:), 

20F ^/2 
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we obtain the following recurrence for mg: 
(4.22) 

1 ^ fs\ T{ka' + i)r((s - k)a' - i) sT{sa' - 1) 

To obtain the form of the recurrence in (|4.2()(l . write H4.22|l in the form 



^ s — 1 ^ s— 1 

-j^ ^ es,fc + = -^-^ '^^{es,k + e.s,s-k) ■ 



2V^ 

^ k=\ " k=l 

and simplify. □ 

Remark 4.3. In going from (|4.22l) to (|4.20() we symmetrized by collecting coefRcients 
of mkrus^k- We might also have symmetrized from the start by choosing the 
splitting probabilities as 

Pn,k := ^{Pn,k +Pn,n-k)- 

In the particular case of Cayley trees this leads to the same splitting probabilities 
as for the Union-Find recurrence studied in (121 1111 H]. 

We can now show convergence in distribution via the method of moments. 

Theorem 4.4. Let a > 1/2. Define ct^ := t^<^"{t)/<^{t) and a' -.^ a + \. Then, 
as n oo, 

with convergence of all moments, where X'"^ has the unique distribution whose sth 
moment rUg = ms{a) is given by 

mi = — = 

V2T{a) 

and for s > 2 by the recurrence H4.20() . 

Proof. One need only check that the m^'s satisfy Carleman's condition. This has 
already been established in |S] . □ 

Remark 4.5. It is curious that a^^n^" X„ has the same limiting distribution as 

an-"' ^ |r„r. 

Here T is a random simply generated tree and \Ty\ denotes the size of the tree 
rooted at a node v. This was established in 0. 

For the case < a < 1/2 it is convenient instead to consider the random variable 

An := An - fJ.n, /i := — j=. 

[Note that, by (|4.13l) . /in is the lead term in the asymptotics of EX„ when a < 1/2.] 
Using 

(4.23) Xn^XK„+X:_K^+tn, u > 2; X^^ti-^i. 
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Define fin ■— EX^ and 



n>l 



Then, in analogous fashion, H3.5|l - H3.10|l hold with /i['*l(z) replaced by p}-'^\z) 
Observe that, by (|4.11|1 . 



(4.24) 



{Z-°'+i + + 1 + Z^/^)A. 



We can use H4.24|l and (|3.5|) - H3.10|) to show that Proposition ITT] holds for a < 1/2 
with /il^l (z) replaced by /xl"! (z) and q changed to 2a — e, for sufficiently small e > 0. 
It follows then that Xn — ^in has (after scaling) a limiting distribution. 

Theorem 4.6. Let a < 1/2. Define cr^ r2$"(T)/$(r) anrf a' := a + i. Then, 
as n oo, 

with convergence of all moments, where X^"^ has the unique distribution whose sth 
moment rus = ms(oi) is given for s = 1 by 

r(a-^) 

mi — — ■= — 

V2r(a) 

and for s > 2 by the recurrence 1)4.20(1 . 

Finally we turn our attention to the case a = 1/2. Now, we define 

Xn Xn -j^nlnn — Sn with d ^ A. H i^il + 21n2), 

\/2tt cJtt J2tt 



with Li defined at (HHj). Then [cf. llTT|l ] 



-'^n — A:„ + -'^n-AT,, + J ^ > 2, 



with Xi = 1 — (5 and 



tn.k • — 



/27r 



fc In A; + (n — fc) In (n — fc) — n In n + 



2tt 



-n 



1/2 



cr 



As in the case a < 1/2, it is easily checked that l|3.5|l - ((3.1()|l hold with /i[*l(z) 
replaced by fJ-^^z) and r|f ' at (|3.4|) replaced by 

n-l 



(4.25) 



E 



•S1+S2 + S3=S 
S2,S3<S 



Si, S2, S3 



-k- 



k=l 



The limiting distribution is given by the following result. 
Theorem 4.7. As n oo, 



where tuq = 1, mi = 0, and for s > 2, 



r(s-i) 



20Fr(i 



E 



2' Si+S2+S3=S 



Si, S2, S3 



/27r 



rris'i'^Si J Si, S2, S3 J 
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with 

Js^.s^.s^:^ / [xlTix+{l-x)liiil-x)Y'x''^-^{l-xy^-^ dx. 
Jo 

Consequently 

where X'^^^^^ has the unique distribution whose sth moment is given by nig. 

Proof sketch. We provide an outline of the proof, leaving the details to the reader. 
We claim that it is sufficient to show that 

(4.26) p"/tlr„ = [C.+o(f)]n^~i, 

with Co = c, Ci = 0, and for s > 2, 

Indeed, defining ms '■= o'^'^c^ Cs and proceeding as in Theorem 1321 yields the 
claim. 

To show (|4.26|) . we proceed by induction. The case s = is easily checked, and 
the case s = 1 follows from (|4.17|) . For s > 2 we use the induction hypothesis and 
approximation of sums by Riemann integrals in H4.25|l to get 

where 

^ I I I C ' , 



^ \ Si,S2,S3/ Vv27r ' 



S2,S3<S 

Since we know a priori that f I^I (z) is amenable to singularity analysis it follows that 
[cf. (ICTl ] 

and completing the computations as in the proof of ProDOsition l4.1l vields the proof 

of Esnii- □ 

5. One-sided destruction 

5.1. Expectation. We study equation H3.17|l for the toll t„ — n" with a > 
and start by establishing a singular expansion for the expectation Since 
^[ol(z) = T{z), we have from (|3.14|l that 

rill (z) ^ t{z) [aizT'{z)T{z) + aoT^{z)] , 

which has already been considered in Section 14.11 In the remaining part of Sec- 
tion l5.ll we suppose now a ^ |, . . .} U {0, 1,2,.. .}. (The complementary cases 
are covered in the proof of Theorem 15.11 ) Then a compatible singular expansion 
for A^^{z) is available at (|4.9|) . This leads to the expansion (|4.10() for g^^\z) and 
consequently, using (|2.7() . to 

gi^Hz) r(«+i)^__, 



T(z) 2p0F 



.^-a-i ^ z-^'^A + Z-^A + Z-'^^^A + A. 
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Integrating the last expression gives the singular expansion 

r ^ r(a+j) ^,„ ^ z-"+'^A + Z-^+^A + A + Z^'^A. 

Jo T(t) 2a^/TT 

Now using H3.17|l . we obtain easily the desired expansion for 

(5.1) ~ ^r(a + |) ^_„ ^ + + + ^i/^A 

Via singularity analysis, we thus get the following expansion for the coefRcients: 

iVTTl (a + 1) 

which together with (|2.2() yields the full asymptotic expansion 

(5.2) mL'' ~ ^£i^±iL^"+i + „"AA + n"-^AA + AA, 

V2r(a + 1) 

ith cr defined at H2.3(l . 



wi 



5.2. Higher moments and limiting distributions. We state the main result of 
this section: 

Theorem 5.1. Let a > 0. Define a :— '^\J \(r) ' '^"■'^ ^' + |- Then, for 

[si 

toll function tn — n'^ , the moments fin ■— satisfy the following asymptotic 

expansion as n — > cx); 

2^/2llr(,v + i)" ^' 

_fi-e if a £{0,1/2} 

|min{Q!, 1/2} otherwise, 

where e > is arbitrarily small. Thus the normalized random variable Yn converges 
weakly to a random variable y (") ; 

where y(") has the unique distribution with (for s > 1) sth moment 



Os/2 11 



r(ja') 



2V2lir(ja' + i)- 

/n particular when a = (i.e., tn = Ij, CTn^^/^y„ converges weakly to a standard 
Rayleigh distributed random variable Y'^^^ with density 

f{y) = ye-y"^^, y > 0. 
[si 

In this case the asymptotics of fin can be sharpened to 

' log n ^ 



l^n 



2s/2r(s±l) 



n 



l + O 



18 



JAMES ALLEN FILL, NEVIN KAPUR, AND ALOIS PANHOLZER 



Proof. We use induction on s. We begin with a > 0. Observe that it is sufficient to 
show that the generating functions /i[''l(z) admit the asymptotic expansions (|5.3|l 
around their dominant singularities at z = p. Then, using singularity analysis, the 
claim follows. What we will show is that 

(5.3) m[^1(z) = n ^ 7.^'^°'+Ho(|Z|— '+H^), 

where 

n;=ir(j«o 

First we consider s = 1, where we immediately obtain from the full expan- 
sion H5.1|l that 15.3|l is true for all a ^ {i, |, . . .} U {1, 2, . . .}. If on the other hand 
a £ {^, |, . . .} U {1, 2, . . .}, then, repeating the computations of Section [Ql it is 
easily seen that logarithmic terms appear in the expansion of ii^^^z). But apart 
from the case a = ^, they don't have an influence on the main term or on the 
asymptotic growth order of the second-order term. If a = i , one observes that the 
general formula for the main term holds, but the bound for the remainder term 
is different: 0(|logZ^^|), not 0(1). Summarizing these cases, the expansion (|5.3(l 
holds for s = 1. 

Next we assume that H5.3|l holds for all 1 < S2 < s with a given s > 1. From H5.3|l 
follows the expansion 

2 2 y/rrp 

which holds for all 1 < S2 < s. Together with fi^'^^z) = T{z) and aiT = 1, this 
gives the following singular expansion: 

T{z){aiZ^,^i^'-Hz) + aofi^'^^z)) = 

Under the assumptions si + S2 ~ s and S2 < s, we get via singularity analysis 
the expansion 



[nz){a,zd,f,^^^\z) + ao^^^'^\z))] 



[Uj2<=2 + i)/2y^Flli=l r(ja' + i)" +U{n ), i S S2 < s. 

Thus under the assumptions given above, the dominant contribution to r['*l(z) 
is obtained when S2 = s — 1 and si = 1, giving the expansion 

P 2«/2V^ n r(ja' 4- i) + ^ ^' 

which in turn yields the following singular expansion for r['*l(z): 
(5.4) rW(z) = i^lJ^^.Z--' +0(|zr-'+^). 
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Immediately from (|3.16|l and H2.6|l follow the expansions 

9 ^''> 2(-+i)/2/^p^'' +^u^i ; 



and 

T{z) ~ W+^^p 
Integrating leads to 



/o T{t) 2(«+i)/20FGsa' - i; 

Using (|3.17() and H2.1() . we obtain H5.3() and Theorem 15. II is proved for a > 0. 

The case a — Q has already been proved in 15 , where the distribution has been 
characterized by its moments. Therefore we describe only very briefly how to obtain 
this result with the present approach. 

One need only show by induction the singular behavior 

(5.5) m'^I(^) 

ff lnZ-i+yl + 0(|Zi/2logZ-i|), s = l, 

= i V2arZ-i/2+o(|iogZ-i|2), .s - 2, 

[^a^"i2(^-i)/2r(f + i)r(^)z-(^-i)/2 + o(|z-t+iiogZ-i|), s > 3. 

The desired result then follows by applying singularity analysis and the duplication 
formula for the F-function. 

To begin the proof of (|5.5() . first we remark that for s — 1 one proceeds as in 
Section ism and gets the full expansion 

{z) = ^\nZ-^+A+ (Zi/2 log Z-^)A + Z^I^A + (Z log Z-^)A, 

which of course gives ()5.5|l in that case. Assuming that H5.5|l holds for 1 < S2 < s 
with a given s > 2, we have the singular expansion 



d,yy'-\z) - ^a^2-i2(^2-i)/2r f fa + 1) r 

\jTxp V 2 / 



Si + 1 



2-(s2 + l)/2 



2 

+ 0(|Z-^logZ-i|). 

Under the restrictions s\ + S2 = s and si < s, we obtain via singularity analysis 
the expansions 

I ( 'Ji®'^^) [T(z)(aiz9,A^[^^l(z) + aofi^'^K^))] } 

^ _^^3.~i2^r (I + l) + 0(n*~i logn) 

for 1 < S2 < s. [For S2 — 0, an expansion is already available at 1)4.7(1 .] Under the 
given restrictions, the dominant contribution to r^'*! (z) is obtained when S2 = s — 1 
and we obtain the following singular behavior of r['*l(z): 



rW(-) = ^-^-^2t-iF (^i±i) F (f + 0{\Z-^ \ogZ-'\). 
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We proceed with 



and, due to (|3.17|) and l|2.1|l . integrating gives H5.5(l for s > 2 and completes the 
proof. □ 
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